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The B° lifetime has been measured with a sample of 23 million BB pairs collected by the BABAR 
detector at the PEP-II e + e~ storage ring during 1999 and 2000. Events from the semileptonic decay 
B° — > D*~£ + ve have been selected with a partial reconstruction method in which only the charged 
lepton and the slow it from the D*~ — > D°ir~ decay are reconstructed. The result is 

t b o = 1.529 ± 0.012 (stat) ± 0.029 (syst) ps. 
PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The technique of partial reconstruction of D*~ mesons 
(charge conjugate states are always implied), in which 
only the slow pion from the D*~ — > D°ir~ decay is re- 
quired, has been widely used in the past Q to select large 



samples of reconstructed B mesons. This technique pro- 
vides a way to measure the combination of CKM angles 
(2/3 + 7) w itb B° -> D*-tt+ decays f|. The application 
of this method to the semileptonic decay B° — > D*~£ + h*i 
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allows the method's validation, while providing a tool 
for precise measurements of several properties of the B , 
including its lifetime, t b o , and the B°B° mixing param- 
eter, Am^. A precise measurement of t#o is presented 
herein. 

The data used in this analysis, recorded by the BABAR 
detector at the PEP-II storage ring during 1999-2000, 
correspond to an integrated luminosity of 20.7 fb _1 col- 
lected on the T(AS) resonance (on-peak events) and 2.6 
fb _1 collected 40 MeV below the resonance (off-peak) 
for background studies. Samples of simulated BB events 
were analyzed through the same analysis chain as the real 
data. The equivalent luminosity of the simulated data is 
approximately equal to the on-peak data. 

A detailed description of the BABAR detector and the 
algorithms used for track reconstruction, particle iden- 
tification, and selection of BB events is provided else- 
where ||; a brief summary is given here. Particles with 
momenta p>,170MeV/c are reconstructed by matching 
hits in the silicon vertex tracker (SVT) with track ele- 
ments in the drift chamber (DCH). Since lower momen- 
tum tracks do not leave signals on many wires in the 
DCH due to the bending induced by the magnetic field, 
they are reconstructed by the SVT alone. Electrons are 
identified with the ratio of the track momentum to the as- 
sociated energy deposited in the calorimeter (EMC), the 
transverse profile of the shower, the energy loss in the 
drift chamber, and the information from the Cherenkov 
detector (DIRC). The efficiency for electron identifica- 
tion in the acceptance of the electromagnetic calorimeter 
is about 90%, with a hadron misidentification probability 
equal to 0.15%. Muon candidates are required to have 
a path length and hit distribution in the instrumented 
flux return and energy deposition in the EMC consistent 
with that expected for a minimum-ionizing particle. The 
Cherenkov light emission in the DIRC is then employed 
to further reject kaons misidentified as muons, by requir- 
ing muon candidates to have a kaon hypothesis probabil- 
ity less than 5%. These criteria yield 74% muon efficiency 
with 2.6% hadron misidentification probability. 

Semileptonic B° decays are then selected by search- 
ing for the high momentum charged lepton (£ — e,/i) 
from the B° decay and the slow pion (tt s ) from the 
D*~ — ► D°tt~ decay. To reject leptons from semileptonic 
charm decay and misidentified hadrons, the momentum 
of the lepton candidate in the T(AS) rest frame (p|) is 
required to be in the range 1.4 < p* t < 2.3 GeV/c; that of 
the ir s (p* ) has to be less than 0.19 GeV/c. The kine- 
matics of the decay are exploited for further background 
suppression as follows. As a consequence of the limited 
phase space available in the decay D*~ — ► D°it~ , the ir s 
is emitted within a one-radian wide cone centered about 
the D*~ direction in the T(4S) rest frame. The D*~ four- 
momentum can therefore be computed by approximating 
its direction as that of the 7r s , and parameterizing its mo- 
mentum as a linear function of the n s momentum, with 



TABLE I: Composition of the data sample in the signal re- 
gion. The error is the sum in quadrature of the statistical and 
systematic errors from the fit to the data M v 2 distribution. 



Sample 


# of events 


Fraction (%) 


Signal Region 


172,700 




Backgrounds 






Continuum 


19, 600 ± 400 


11.4 ±0.2 


BB comb. 


52, 700 ± 1,400 


30.0 ±0.8 


B+ 


8, 700 ± 4, 400 


5.0 ±2.5 


B" signal 


91, 700 ±4, 600 


53.6 ±2.6 



parameters obtained from the simulation. The neutrino 
invariant mass can be computed from the four-momenta 
of the B°, D*~ , and I with the relation 

M 2 = {P B o -P D .- -P e ) 2 . 

The momentum of the B° in the T(4S) rest frame, on av- 
erage 0.34 GeV/c, is neglected. Mj 2 peaks approximately 
at zero for signal events, whereas background events are 
spread over a wide range. 

The B° decay point is determined from a vertex fit 
of the 7r s and I tracks, constrained to the beam spot 
position in the plane perpendicular to the beam axis (the 
x-y plane). The beam spot is determined on a run- by-run 
basis using two-prong events || . Rs size in the horizontal 
direction is 120 /xm. Although the beam spot size in 
the vertical (y) direction is only 5.6 /xm, a beam spot 
constraint of 50 /im is applied to account for the flight 
of the B° in the x-y plane. Only events for which the \ 2 
probability of the vertex fit, Vvi is greater than 0.1% are 
retained. 

A selection is applied on the combined likelihood for 
Pi , , and Vy , which results in a signal-to-background 
ratio of about one to one in the signal M v 2 region, defined 
as M 2 > -2 GeV 2 /c 4 . Figure @ shows the M 2 distri- 
bution of data events used to measure Tgo when the I 
and the ir s have opposite-sign charges. Same-sign events 
are used as a background control sample. The individual 
distributions shown in Fig. |l| are obtained by fitting to 
the data the contributions from continuum events, ob- 
tained from the off-peak data, and from BB combinato- 
rial background, B° signal, and B + resonant background, 
as predicted by the simulation. The B + resonant back- 
ground is due to intermediate production of higher mass 
charm resonances (denoted as D**). The fit determines 
the composition of the selected sample, which is reported 
in Table || for the events in the signal region. 

The PEP-II collider produces BB pairs moving along 
the beam axis (z direction) with an average Lorentz boost 
of (/?7) = 0.55. Hence, the two B decay vertices are sep- 
arated on average by (Az) i=s 255 /im. The position of 
the B° — > D*~£ + ve ("decay") vertex is reconstructed as 
described above. The decay point of the other B is de- 
termined from a selection of the remaining tracks in the 
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FIG. 1: The M„ 2 spectrum of the selected events. The data 
are represented by the dots with error bars. The result of the 
fit with shapes from the simulation are overlaid. 



event using the following criteria. In events that have an- 
other lepton with momentum p* > LlGeV/c, the other 
B vertex is computed with only this lepton track con- 
strained to the beam-spot in the x-y plane. Otherwise, 
all the tracks with a center-of-mass angle greater than 
90° with respect to the 7r s direction are considered. This 
requirement is used to remove most of the tracks from the 
decay of the D° daughter of the D*~ , which would oth- 
erwise bias the reconstruction of the other B vertex posi- 
tion. Simulation shows that in about 75% of signal events 
the other vertex has no tracks from the D° decay. The 
selected tracks are then constrained to the beam-spot in 
the x-y plane. The track with the largest contribution 
to the vertex x 2 , if greater than 6, is removed and the 
fit iterated until no track fails this requirement. Vertices 
composed of just one track that is not a high momen- 
tum lepton are rejected in order to reduce the number 
of poorly measured vertices. The lifetime is determined 
by measuring the quantity Az — Zdccay — Mother, where 
Zdecay (bother) is the position along the beam line of the 
decay (other) vertex. The proper time difference is then 
computed with the relation At = Az/(c/3j). A fit with a 
double Gaussian to the At residuals in the Monte Carlo 
simulation shows that one half of the events are contained 
in the narrower Gaussian, which has a width of 0.7 ps. To 
remove badly reconstructed vertices, all events for which 
either \Az\ > 3 mm or ctaz > 500 /im are rejected, where 
a Az is the uncertainty on Az computed for each event. 

Tgo is obtained from a binned maximum likelihood 
fit to the two-dimensional At, a At distribution. To 
save computation time events are grouped in a two- 



dimensional space consisting of 100 At and 25 a At bins. 
The At distribution of signal events, J 7 (At, er^t, Tgo), is 
described by the convolution of the decay probability dis- 
tribution 



/(A* t rue 



i 



2r no 



ea;p(-|Attruc|/TBo), 



with the experimental resolution function, which is 
parametrized by the sum of three Gaussian distributions. 
The two narrow Gaussians, which account for more than 
99% of the events, have the form 



g(S(At),a At ) 



1 



2nS(JAt 



(S(At)-b)\ 



where 6 (At) — At — Attrue is the difference between the 
measured and the true value of At, b is a bias due to 
the charm tracks in the other vertex and resolution ef- 
fects, and the scale factor S is introduced to account for 
possible misestimation of the calculated error a At on the 
proper time difference. The third Gaussian of fixed bias 
(— 2ps) and width (8ps) is added to account for badly 
measured events ("outliers"). 

B + background events that peak in the M„ 2 signal re- 
gion are described by an identical function, with the same 
resolution parameters as for the B° signal events, and an 
effective B + lifetime of 1.57 ps. This value, obtained by 
fitting simulated B + events, is smaller than the value of 
1.655 ps generated in the simulation due to D° tracks 
from the decay vertex being included in the other vertex. 

The At distribution of continuum background events 
is modeled as the sum of two components, one with non- 
zero lifetime and the other with zero lifetime, convolved 
with the same single Gaussian resolution function. The 
parameters of the resolution function, as well as the life- 
time and the fraction of events with non-zero lifetime, 
are all determined with the off-peak events that satisfy 
the selection criteria. 

The At distribution of the combinatorial BB back- 
ground is modeled as the sum of a non-zero and a zero- 
lifetime component, with a resolution function that is 
the sum of three Gaussians. All parameters are deter- 
mined from the data by fitting the measured At distri- 
bution of the events in the sideband region, defined by 
-10 < My 2 < -4 GeV 2 /c 4 . The Monte Carlo simulation 
shows, however, that there are small differences in the 
lifetime and in the fraction of events with non-zero life- 
time between the signal region and the sideband. These 
differences are also observed in the data by separately fit- 
ting signal region and sideband events in the same sign 
£tt s background control sample. The results from the like- 
sign fits are used to scale the two background parameters 
from the sideband to the signal region. 

The function used to fit the data is the weighted sum 
of the four contributions: 
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F(At,a At \T B o) = [1 - f B+ {M v 2 ) - f c {M v 2 ) - f Bn (M u 2 )}F B o(At,a At ,T B o) + f B+ (M„ 2 ) T B+ (At, a At) + 
f c (M u 2 ) T c (At, <TAt) + ! BB {M 2 ) T BB (At, a At), 

I 



where the functions T B o , T B + , J- c and T BB describe 
the measured decay time difference distributions for the 
signal, peaking B + , continuum, and BB combinato- 
rial background, respectively. f B +, f c and f BB are 
the probabilities that the event is from the B + , contin- 
uum, or BB background, computed for each event on 
the basis of the measured value of M y 2 . Simultaneously 
with t b o, the following parameters of the signal reso- 
lution function are fitted: the scale factor of the first 
Gaussian, Si = 1.02 ± 0.02, the scale factor of the sec- 
ond Gaussian, 5 2 = 2.4 ± 0.1, the bias of the first Gaus- 
sian, bi = —0.120 ± 0.009, and the fraction of outliers, 
fo = (0.2 ± 0.1)%. The fraction of events contained in 
the second Gaussian, / 2 , and its bias 62 are fixed to 7% 
and —0.85 ps, respectively. 

The result of the fit is t^ w = 1.482 ± 0.012 ps, where 
the error is statistical only. Figure [2] shows the com- 
parison between the measured At distribution and the 
fit result. The probability of obtaining a lower likeli- 
hood, evaluated with a Monte Carlo technique, is 18%. 
This raw lifetime must be corrected for the bias in- 
duced by the tracks from the D° that are not rejected 
by the n s cone cut. A multiplicative correction factor 
of K Eo = 1.032 ± 0.007(stat) ± 0.007(syst) is computed 
from the simulation. The statistical error arises from 
the number of simulated events. The dominant system- 
atic uncertainty corresponds to the full variation in TZ B a 
(0.66%) obtained by smearing the At resolution in the 
simulation to match that in the data. A second system- 
atic uncertainty is computed by comparing in data and 
simulation the fraction of charged tracks from D° decays 
outside the 7r s cone for a subset of events in which the 
D° is fully reconstructed in the K + ir~ , K + ir~ir°, and 
_ftT + 7r~7r + 7r~ final states. The maximum discrepancy be- 
tween data and simulation corresponds to a variation of 
±0.24% in the value of 1Z B0 - The corrected value of the 
B° lifetime is then 



T B o 



raw'n _ 
~B° '^D° 



1.529 ± 0.012 ps. 



The systematic error on t b o is computed by adding in 
quadrature the contributions from several sources, de- 
scribed below and summarized in Table p[ 

The fractions of B + , continuum, and combinatorial 
BB events are varied by the uncertainties obtained from 
the Mj 2 fit (see Table |). The parameters of the contin- 
uum and combinatorial BB At distributions are varied 
by their uncertainties, accounting for their correlations. 
As described above, the fraction of events with non-zero 
lifetime and the lifetime of the combinatorial BB back- 
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FIG. 2: At distribution for selected events in the data (points) 
in linear (upper) and logarithmic (lower) scale. The lifetime 
fit result is superimposed on the data. The hatched his- 
tograms show the contributions from the background sources 
described in the text. 



ground computed from the sideband are corrected with 
the same-sign control sample. This method is validated 
by a simulation study, and the statistical error of the val- 
idation is included in the background systematic error. 
The effective B + lifetime is varied by ±3%, which is the 
sum in quadrature of the world average error on the B + 
lifetime and the statistical and systematic uncertainties 
of the D° bias correction. 
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The parameters of the signal resolution function that 
are not determined in the fit to the data are varied within 
conservative ranges (/ 2 between 0.03 and 0.13, and b 2 be- 
tween —1.5 and ps). Several different analytical expres- 
sions are used to represent the small fraction of outliers. 
The fit is also performed by allowing the scale factor and 
the bias of the narrow Gaussians to depend linearly on 
a At or on the lepton polar angle. The maximum change 
with respect to the result with fixed parameters is taken 
as the systematic error due to the parametrization of the 
resolution function. 

The bias due to the event selection is found to be com- 
patible with zero by fitting with an exponential function 
the true proper time difference of signal events selected in 
the simulation. The statistical error of this test is added 
to the systematic error. 

The statistical and systematic errors on TZj^o are prop- 
agated to the final error. A possible bias induced by 
the presence of tracks from charm decays produced by 
the other B meson is investigated in the simulation by 
varying within their uncertainties the relative fractions of 
charmless, single charm, and double charm events, and 
also by varying the relative fractions of D + ,D°, D s , and 
A c hadrons. The z length scale is determined with an 
uncertainty of 0.4% from secondary interactions with a 
beam pipe section of known length. The dependence of 
the result on several different variables (angular width 
of the 7r s cone used to reject D° tracks, soft pion mo- 
mentum, lepton momentum, polar and azimuthal angle, 
alignment conditions) is carefully inspected; no statis- 
tically significant effect is found. No difference in the 
result is observed if t b o is determined with an unbinned 
maximum likelihood fit. A final relative error of ±1.9% 
is found by adding in quadrature the uncertainties from 
the above sources, as listed in Table ||. 

In conclusion, a sample of about 92000 B° — > D*~£ + U£ 
decays is selected by partial reconstruction of the 
D*~ — > D°tt~ decay. It is used for a measurement of 
the B° lifetime. The value obtained, 

t b o = 1.529 ± 0.012 (stat) ± 0.029 (syst) ps, 

is consistent with a recent BABAR measurement Qj and 
with the world average || . It is currently the most precise 
single measurement of this quantity. 
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TABLE II: Contributions to the systematic error. 



Source 


0V rO /t s o(%) 


Continuum fraction & parametrization 


0.36 


BB fraction & parametrization 


0.68 


B + fraction & parametrization 


0.64 


Resolution model 


1.14 


Event selection bias 


0.30 


D° bias {K B o) 


0.95 


Bias due to charm from the other B 


0.21 


z scale 


0.40 


Total 


1.89 



